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ABSTRACT 

We extend a recently developed galaxy morphology classification method, Quantita- 
tive Multiwavelength Morphology (QMM), to connect galaxy morphologies to their 
underlying physical properties. The traditional classification of galaxies approaches 
the problem separately through either morphological classification or, in more recent 
times, through analysis of physical properties. A combined approach has significant 
potential in producing a consistent and accurate classification scheme as well as shed- 
ding light on the origin and evolution of galaxy morphology. Here we present an anal- 
ysis of a volume limited sample of 31703 galaxies from the fourth data release of the 
Sloan Digital Sky Survey. We use an image analysis method called Pixel-z to extract 
the underlying physical properties of the galaxies, which is then quantified using the 
concentration, asymmetry and dumpiness (CAS) parameters. The galaxies also have 
their multiwavelength morphologies quantified using QMM, and these results are then 
related to the distributed physical properties through a regression analysis. We show 
that this method can be used to relate the spatial distribution of physical properties 
with the morphological properties of galaxies. 

Key words: galaxies: classification - galaxies: formation - galaxies: evolution. 



1 INTRODUCTION 



The Hubble tuning fork (Hubble 19261 is one of the first 



and most widely used galaxy classification schemes. Even 
though it has provided many insights into the evolution of 
galaxies, their morphologies and other physical properties, 



van den Bergh 1998 I it remains subjective, requiring ex- 



perts to manually classify galaxies, and is not directly re- 
lated to the physical properties of the galaxies. The use of 
one wavelength («450nm) also restricts the Hubble scheme 
and studies such as Block & Wainscoat ( 1991 1 and Jarrett 



(20001 show clear examples of why a multiwavelength ap- 



proach is vital for robust and thorough classification. Many 
galaxies have been discovered that simply do not fit into the 



Hubble scheme (van den Bergh 1976 Sandage & Brucato 



1979). This is further exemplified at high redshift, where 



galaxies are at earlier stages of their evolution. The inclina- 
tion of galaxies also plays a vital part in the classification 
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process, especially for those galaxies that fall into the spiral 
sequence. 

Galaxy classifications need to be compact and physi- 
cally motivated and the above limitations have motivated 
a multitude of classification schemes to achieve these goals 
while addressing the various drawbacks of the Hubble tun- 
ing fork. De Vaucouleurs (1959) extended Hubble's scheme 
by introducing more divisions within classes. For instance, 
spirals were divided into more refined classes based on the 
presence of bars and rings around the galaxy and the spi- 
ral arms were divided into three new classes. The intrinsic 
basis of the Hubble system remains, though, limiting the 
utility of such schemes in progressing towards a connection 
between morphology and the underlying physical processes. 



Simpler classification systems by Morgan ( 1958 1959 1 us 



ing the central concentration of light and a system based 
solely on the stellar population of galaxies by Morgan & Os- 
terbrock ( 1968 I have been successful at providing high dis- 



crimination of galaxies between classes. Kormendy fc Bender| 
( |1996[ ) aimed to revise the original classification of ellipticals 
by Hubble as this is mainly correlated to the inclination of 
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galaxies rather than any intrinsic properties of the galaxies. 
These authors developed two main classes for ellipticals and 
several branches within each type. As a subjective scheme 
however, the classes are not robust and many galaxies have 
properties that belong to multiple classes or none at all. |Ko-| 



ideal candidate to be used in the morphological classification 



rmendy & Bender ( 1996 1 attribute these inconsistencies to 
heterogeneous formation histories. 

More recently automated classification schemes such 



as Artificial Neural Networks (ANNs) (Nairn et al. 1997a 



Bazell 2000 Odewahn et al. 2002 1 have been used to ac- 



commodate the vast quantities of galaxies that require clas- 
sification and to a large extent ANNs eliminate any subjec- 
tive bias in the classification process, as well as being very 



accurate (Ball et al. 2004). ANNs cannot create a new classi- 
fication system but rather replicate visual classification and 
with much higher consistency. These systems require a "test 
sample" previously classified by a human expert on which to 
base classifications which has the disadvantage of allowing 
the same human biases and flaws of the classification system 
to propagate. The use of Self Organising Maps by Nairn et 
al. ( 1997b I however, eliminates the need for a test sample 



and any human influence in the classification process. 

Photometric decomposition techniques which analyse 



the observed distributions of photometric intensity (Simard 
at al.|2002 Peng et al.|2002[ ) and Fourier analysis techniques 



that quantify luminosity distributions of galaxies ( Odewahn 



et al. 2002 Trinidad 1998k have had moderate success in 
II - I I II r 

differentiating galaxies into their respective classes. 

Non-parametric approaches such as the CAS classifi- 
cation scheme have had success in objectively separating 
galaxies into Hubble's classes as well as being applicable at 
high redshifts (z=3) ( |Conselice|2003| | Abraham et al.|1996 |. 
The scheme uses 3 properties, concentration, asymmetry 
and dumpiness, which quantify aspects of galaxy morphol- 
ogy. These quantities identify formation histories, merging 
activity and areas of high star formation activity ( |Conselice| 
2003). This technique can be easily applied to the decom- 



position of the distribution of other physical properties in 
galaxies. The inclusion of additional parameters such as a 
Gini Coefficient has been shown to produce more refined 
separations of galaxies but at the expense of increasing the 
dimensionality of the classification scheme ( |Lotz et al.|2004[ 
2006). 



Classification of galaxies by physical properties alone 
has not been extensively carried out, although new tech- 
niques such as Pixel-z have emerged that enable the extrac- 
tion of information about the physical properties of galaxies 
by fitting spectral energy distribution (SED) templates de- 



rived from stellar population evolution ( |Bruzual fc Chariot 
|2002||Conti et al. 12003] ). Each pixel is fitted with templates, 
giving a localised analysis of the physical properties within 



the galaxy (Welikala et al.||2007 2008 2009) 



Shapelet decomposition promises a new approach in 
the morphological classification of galaxies. Shapelets are 
Gaussian- weighted Hermite polynomials (Rcfrcgicr 2003). 
They are also the eigenstates of the Quantum Harmonic Os- 
cillator (QHO) Hamiltonian, and are thus well understood 
( jRefre gicr 2003). They have been shown to be useful in im- 



mcasurcmcnts 



age simulation ( Masscy et al. 2004 1 and gravitational lensing 
( Chang & Refregier 2002 Rcfrcgier & Bacon 



|2003[ ). Shapelets use all the information about the shape 
of a galaxy and form a complete set thus making them an 



of galaxies ( Kelly & McKay 2004 1 . Shapelets are a central 



component in a new objectively developed and automated 
classification system known as the Quantitative Multiwave- 
length Morphology (QMM). QMM uses shapelets to decom- 
pose the galaxy images and a Principal Component Analysis 
(PCA) to reduce the dimensionality of the data followed by 
a Mixture-of-Gaussian models to objectively identify partic- 
ular morphological classes of galaxies ( Kelly &: McKay|2004 



2005 1 . Shapelets are not a compact form of classification and 
require a PCA to account for this. The technique uses galaxy 
images in multiple filters, and currently images in the Sloan 
Digital Sky Survey (SDSS) filter set (ugriz) has been used. 
Kelly fc McKay| ( |2004[ ) and |Kelly fc McKay| ( [2005] ) show that 
this technique consistently reveals previously established re- 
lationships such as that between Hubble type and colour, 
as well as broad connections between morphology and the 
physical properties of galaxies. 

We aim to identify relationships between the physical 
properties, measured with Pixel-z and quantified using CAS, 
and the morphological properties quantified using QMM. 



2 QUANTIFICATION OF PHYSICAL 
PROPERTIES 

Our initial objective is to extract the physical properties 
from galaxy images and then quantify the spatial distribu- 
tions of these physical properties. The tools used for this 



process are Pixel-z (see Welikala et al. 2007 2008 20091, 



for extracting the physical properties, and CAS (Conselice 
|2003[ ), to quantify their spatial distribution (the collective 
process will be known as Pixel-z : CAS herein). 

These quantities are then compared with the results 
of the QMM analysis of the same galaxy images, through 
a regression analysis, to analyse how well QMM describes 
the spatial distribution of the physical properties in galax- 
ies. This indicates the extent to which we can use QMM to 
connect physical and morphological properties of galaxies, 
providing us with the possibility of developing a compre- 
hensive galaxy classification scheme that incorporates both 
physical and morphological properties of galaxies. 

The data was obtained from the fourth data release 
(DR4) of the Sloan Digital Sky Survey (SDSS). The SDSS 
used a dedicated 2.5m telescope located at Apache Point 
Observatory in New Mexico, USA together with a 142 
megapixel camera in drift-scan mode to obtain images and 
spectroscopy over about a quarter of the sky (for more de- 



tails, 


see York et al 2000| |Hogg et al. 


2001; |Smith et al. 


2002 


ISubbaRao et al. 20021 IPier et al. 


|2003| |Ivezic et al. 


2004 


|Tucker et al.|2006 


) 



2.1 Parameterization of Pixel-z Output 



Using the Pixel-z output from Weli kala et al.| ( |2008[ ) we 
quantified the distribution of the physical properties in 
galaxies through the CAS technique developed by |Conselice| 
(2003). The CAS analysis uses the parameters concentra- 



tion, asymmetry and dumpiness (C, A and S, Conselice 



20031. We calculated these three parameters for four phys- 
ical properties, age, star formation rate, colour excess and 
metallicity, each of which were obtained from the Pixel-z 
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analysis of Welikala et al. (2008). The same process was 



carried out for the r-band (616. 5nm) image of the galaxy. 

From the r-band CAS parameters we can contrast the 
relationships between QMM and physical properties with 
that between QMM and the photometric distribution for 
the galaxies. The r-band is convenient for this process as it 
lies in the middle of the SDSS filters and typically shows a 
high signal to noise ratio. 

The concentration parameter is calculated from the ra- 
tio of the radii containing 20% and 80% of the total flux. 
Asymmetry is derived from subtracting the image of a 
galaxy rotated by 180 ° from its original image. The dumpi- 
ness is quantified through the ratio of the flux in high spatial 
frequency structures within a galaxy to its total flux. These 
CAS parameters provide a morphological description of the 
physical properties of galaxies by quantifying their distribu- 
tions derived from Pixel-z. 



3 QUANTITATIVE MULTIWAVELENGTH 
MORPHOLOGY 

QMM is a morphological classification method developed by 
|Kelly fc McKay|p054l[2005t , using galaxy images observed 
through the five SDSS filters (u, g, r, i, z). Our implemen- 
tation of QMM involves two steps, shapelet decomposition 
and a Principal Components Analysis (PC A). The shapelet 
decomposition breaks down the image of the galaxy repre- 
sented in all 5 filters into a set of mathematical functions 
which are quantified by a set of coefficients. To adequately 
reconstruct galaxy images the number of coefficients can be 
in the hundreds, so PCA is used to reduce the number of 
coefficients to a more manageable size. The outputs of the 
QMM analysis are the final PCA components, which encom- 
pass the majority of the variance in the data. 

For the shapelet decomposition 37876 galaxy images 
from the Pixel-z:CAS analysis were available, comprised of 



the volume-limited sample analyzed by Welikala et al. ( 2008 
2009). The resulting data set excluded 432 galaxies, the im- 
ages for which contained numerous zero flux values that 
could not be processed in the PCA analysis. The sample 
was further restricted with an 11a cut on the variances in 
each dimension in order to ensure that the PCA was not 
compromised by outlying data points. The final sample for 
the QMM analysis consisted of 31703 galaxies with images 
from all five filters and common to both the Pixel-z : CAS 
and QMM analysis. 

The galaxies in the sample span a redshift range of 
0.00278 < z < 0.231. The analysis was first carried out 
over the entire redshift range. The analysis was then inde- 
pendently applied to four redshift bin sub-samples contain- 
ing approximately equal numbers of galaxies. This approach 
allows us to identify possible biases due to the coarser spa- 
tial sampling of galaxies at progressively higher redshift, im- 
posed by a fixed observational pixel scale, and identifying 
whether it is necessary to artificially redshift the galaxy im- 
ages to mimic a common spatial sampling. 

By comparing the results from QMM with those of 
Pixel-z : CAS we can explore how QMM enables us to con- 
nect galaxy morphology with the underlying physical prop- 
erties. The colour (or flux ratio) information encoded in the 
multiwavelength images used for this technique provides the 



key, as galaxy colours are a consequence of the combination 
of stellar evolutionary processes and multiple stellar popula- 
tions. This suggests that there should be a direct connection 
between a quantitative morphology derived from multiwave- 
length images and the underlying properties of the stellar 
population within the galaxies. Therefore, we aim to deter- 
mine the extent to which QMM is a reliable descriptor of 
the underlying physical properties of galaxies by comparing 
it with the results of the Pixel-z : CAS analysis. 



4 RESULTS 
4.1 Pixel-z: CAS 

The CAS analysis provides morphology descriptions for the 
four sets of physical parameters from Pixel-z: age, star for- 
mation rate, colour excess and metallicity. For each of these 
the concentration, asymmetry and dumpiness parameters 
were measured. 

Our results for the distribution of the star formation 
rate of galaxies in CAS space (Figure [I]) are not consistent 
with those found by Conselice (2003) for the photometric 
images. This is not unexpected, as this is the first time that 
CAS has been applied to a distribution of physical rather 
than photometric properties. The ranges spanned by our 
CAS values are different, extending over a much larger range 



than those of Conselice ( 2003 1. The dumpiness parameter in 



particular, (Figures lb, lc), shows a more restricted distri- 
bution than in Conselice ( 2003 1 , with a number of systems 



having dumpiness very close to zero. The asymmetry pa- 
rameter shows a maximum value of 2, a consequence of the 
mathematical definition of asymmetry. 

Figure [I] (a) shows that at lower concentrations there is 
no relationship between concentration and asymmetry. As 
we move to higher concentrations the asymmetries fall be- 
tween 1.5 and the highest possible value of 2. At high asym- 
metries the star formation activity is occurring only in highly 
localised regions of the galaxy. The proportion of light em- 
anating from these regions within the galaxy compared to 
the central bulge is smaller than in a galaxy with a low 
asymmetry. This leads to a high concentration index. 

dumpiness is not correlated with concentration for 
SFR at any value. Figure [l] (b) shows uniform distribution 
above Concentration = 3 and dumpiness = 0. A possible ex- 
planation for this could be that the dumpiness parameter 
takes into account the galaxy centers as well as the regions 
outside the center. The dumpiness of the galactic center 
would be expected to correlate positively with the concen- 
tration. The dumpiness outside the central regions would 
be expected to correlate negatively with the concentration, 
as the more clumpy the outer regions, the more localised 
the SFR activity. This leads to a higher concentration index 
as in the case of Figure [l] (a). When the dumpiness of both 
these regions are taken into consideration, we are more likely 
to see a large scatter. The lack of correlation between asym- 
metry and dumpiness seen in Figure [I] (c), has a similar 
explanation to that of Figure [l] (b). 

The CAS distributions for other Pixel-z physical prop- 
erties also differ from those for the r-band light distribution 



of galaxies from Conselice ( 2003 1. Simply by analysing these 



relationships in CAS space it is apparent that the distribu- 
tion of physical parameters in galaxies does not mimic its 
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Figure 1. Relationships between the three CAS parameters for the spatial distribution of star formation rate as inferred by Pixel-z, 
with (a) concentration against asymmetry, (b) concentration against dumpiness and (c) dumpiness against asymmetry. 



photometric morphology in a trivial fashion. This highlights 
the fact that a simplistic approach to relating the single- filter 
photometric morphology (parameterised by CAS) with the 
distribution of physical properties within galaxies does not 
show any obvious connection. This in turn supports our mo- 
tivation for using multiwavelength photometric information, 
through QMM, to relate galaxy morphology to the underly- 
ing physical properties. 

The use of CAS in this investigation is not as a morpho- 
logical indicator directly, but instead as a convenient tech- 
nique for parameterising the distributions of Pixel-z physical 
parameters, to allow a comparison with the QMM results 
from the shapelet decomposition. 

It is possible to use the differences in the physical pa- 
rameters to refine the boundaries for the different galaxy 
classes set by Conselice (20031. For instance, Conselice 
( [2003] ) set the CAS ranges for ellipticals to be, concentra- 
tion in the range 4.4±0.3, asymmetry in the range 0.02±0.02 
and dumpiness in the range 0.00±0.04. As we use larger 
numbers of galaxies and at higher redshifts, these values 
will undoubtedly vary. If we use the four physical properties 



and their CAS boundaries for different galaxy types, it may 
be possible to set a more consistent limit for each class of 
galaxies. This extension of CAS to morphology classification 
based on physical properties could be a productive area for 
further investigation, but is beyond the scope of this paper. 

Figure [2] shows the relationship between the physical 
properties with regards to the asymmetry parameter. They 
allow for an understanding of how these properties affect 
each other and how they can be used to define classification 
criteria. Figure [5] (e) shows the best agreement out of all the 
tested parameters as well as the strongest correlation. Both 
colour excess and metallicity agree with the SFR, but the 
correlation is weaker than that between colour excess and 
metallicity. The asymmetry in age does not match with the 
asymmetry in any of the other parameters but the tight- 
ness in the relationship between age and metallicity (Fig- 
ure[2](f)) has potential in being exploited to make finer cuts 
in the morphological classification. Strong relationships be- 
tween these properties allow for establishing criteria for bet- 
ter galaxy classification than if used on their own. 

The spread of galaxies at high asymmetries in Figure [2] 
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is highly constrained. High asymmetries indicate that the 
distribution of the various physical properties will be very 
"patchy". As a result, these properties will be found in 
highly localised regions of each galaxy. As these properties 
have some relation to one another, it is likely that they will 
all be found in the same localised regions, resulting in highly 
asymmetric galaxies having high asymmetries for the mea- 
sured attributes. This also explains the large distributions 
of galaxies at low asymmetries as the distribution of physi- 
cal properties within the galaxy will not be "patchy" allow- 
ing for the bigger spread in the physical properties across 
the galaxy. This spread will be largely independent of the 
other physical properties leading to lower correlations at low 
asymmetries. 

There is a limit of 2 for all the asymmetry parameters 
in Figure [2] This is an artifact of the definition of asym- 
metry. Asymmetry is defined as the absolute value of the 
difference in flux between the original image and the same 
image rotated by 180°, divided by the flux of the original 
image. The maximum possible value for asymmetry that can 
be attained with this definition of asymmetry is two. 

4.2 QMM - Principal Components Analysis 

Performing the shapelet decomposition on the galaxy images 
using a shapelet order of 12 resulted in 91 coefficients per 
filter, making 455 coefficients per galaxy. A lower number of 
coefficients is insufficient to reproduce galaxy images reliably 
and even though a higher number of coefficients would have 
given a more accurate representation of the galaxy images it 
would have been too large to have been effectively reduced 
by the PCA. 

A PCA was used to reduce the dimensionality of the re- 
sults of the shapelet decomposition ( |Karhunen| 1947| |Loeve| 
1978). To carry out the PCA we first used the same pro- 



galaxies, early and late types. Kelly & McKay (20051 de- 



cedure described in Kelly & McKay ( 2005 1 to calculate the 



shapelet coefficients, with two notable differences. As men- 
tioned above, we do not artificially re-sample the images to 
the same redshift. Also, we use a Singular Value Decomposi- 
tion (SVD) technique similar to the one described in |Berry et| 
al. ( 2004 1 to decompose the shapelet coefficients. The SVD 



accounts for the loss of orthogonality in the shapelets that 
result from the pixelization of the images. However, unlike 



Berry et al. (20041, the number of non-zero singular values 



is chosen in an objective manner. In our work, we chose 
the number of non-zero singular values to minimize the esti- 
mated squared error between the true galaxy image and the 
galaxy image reconstructed from the shapelet coefficients. 
We use Stein's Unbiased Risk Estimation to estimate this 



error (Stein| |1981 1, which is easily calculated as a function 
of the number of non-zero singular values. 

The first component contains the majority of the vari- 
ance («75%) in the data set. Even though the subsequent 
components contain a much lower fraction of the variance 
they are nevertheless important. The number of components 
was limited to 8, which contains 92.3% of the total variance 
and is sufficient for our purposes. 

Discrimination between galaxy populations is not ap- 
parent in the first two panels of Figure 3, but becomes clear 
in the third panel, plotting PCA 1 against PCA 4. This bi- 
modality in galaxy morphology has long been known to exist 
and it is related to the two broad morphological classes of 



tected this same bimodality, between their PCA 1 against 
PCA 2. The fact that this bimodality is seen here in a higher 
order PCA component could be a result of many factors. 



Kelly & McKay ( 2005 1 carried out the principal component 



analysis using a sum-of-squares matrix while we used a co- 
variance matrix. Also, even though we used a sigma cut to 
limit any extreme outliers, there may still be fewer distant 
outliers that can significantly affect the results of PCA. Thus 
the location of the bimodality in PCA space does not reveal 
any fundamental property of galaxies as PCA is simply an- 
other way of representing the distribution of galaxies. The 
other separations in PCA space are very similar to the sep- 



arations obtained by Kelly & McKay ( 2005 1 



Similarly, the galaxy types (or combination of them) 
that each PCA component represents will also depend on 
the above factors. To objectively identify particular classes 
contained within the PCAs, for a particular sample, one 
must apply a Mixture-of-Gaussian methods ( |Kelly fe McKay| 



2004 2005) but this falls outside the aims of this experiment. 



This analysis was reproduced for the sub-samples split 
into redshift bins. The distributions of galaxies in PCA space 
in each of the redshift bins are similar to the analysis with 
the full sample. Furthermore, these distributions are consis- 
tent with each other (Figure |4|, as well as those of Kelly & 



McKay ( 2005 1. These distributions of galaxies in PCA space 



in each of the redshift bins fall within very similar ranges. 
These similarities among the different redshift bins allow us 
to conclude that our results are not strongly dependent on 
the treatment of the sample as a whole despite the range 
of redshift. Thus there is no need to simulate a common 
redshift, through artificially redshifting the galaxy images 
by reducing their resolution, as implemented by |Kelly fe| 
|McKay| ( |2005[ ). The principal components for the redshift 
bin subsamples are also correlated against the Pixel-z:CAS 
measurements below, to identify any potential statistical dif- 
ferences. 

As with the Pixel-z : CAS analysis, our objective is not 
to define a classification system through the distribution of 
the QMM parameters of galaxies. Rather we aim to iden- 
tify relationships between the two methods to identify the 
connection between the multiwavelength morphology and 
the spatial distribution of the underlying physical processes. 
The next step in our analysis is to quantify the relationships 
between the two methods. 



5 CORRELATION ANALYSIS 

Morphological classification alone provides a limited ap- 
proach for understanding the properties and evolution of 
galaxies. Using just the physical properties may be more 
straightforward but these schemes do not address the ori- 
gin or evolution of galaxy morphology. However, the multi- 
wavelength nature of the recently developed QMM classifica- 
tion method provides an opportunity to connect morpholog- 
ical classification with underlying physical parameters. The 
colour information encoded in the multiwavelength images 
provides the key, as galaxy colours are a consequence of the 
combination of stellar evolutionary processes and multiple 
stellar populations. This suggests there should be a direct 
connection between a quantitative morphology derived from 
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Figure 2. Comparison of all physical properties against each other in the asymmetry parameter. Distributions are in the same ranges 
for ease of comparison and the line crossing each plot is the one to one relation. Figures (a), (b) and (c) compare the asymmetry of the 
SFR against the asymmetry of age, colour excess and metallicity respectively. Figures (d) and (e) plot the asymmetry of colour excess 
with the asymmetry of age and metallicity while figure (f) compares the asymmetry of metallicity with age. The asymmetry of SFR, 
colour excess and metallicity follow each other to a large extent. 
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Figure 3. Graphs for the first principal component against the next three. Note the different divisions shown. Components 2 and 3 
have a diverging pattern while the 4th component has two clearly distinct groups. These results reproduce well those of Kelly & McKay 
(2005), but the separation for PCA 4 shown above was found for PCA 2 by |KeIIy fc McKay | {2005] . 



multiwavelength images and the underlying properties of the 
stellar population within the galaxies. 

We have used Pixel-z and CAS as a way of quantify- 
ing several physical properties of galaxies so that they can 
be correlated with the QMM method, in order to identify 
how well QMM can represent the distribution of physical 
properties in galaxies. In order to connect the two methods, 
we carry out a regression analysis by performing a multiple 
linear regression fit to identify the extent to which QMM 
correlates with spatial distributions of physical properties 
in galaxies represented through the Pixel-z : CAS method. 
In an ideal case we would find strong correlations between 
the tested parameters, but due to the large scatter seen in 
the above figures this is unlikely. For this reason we focus on 
the relative correlations between the CAS parameters mea- 
sured for the physical properties and the r-band to observe 
whether any physical properties consistently show correla- 
tion coefficients higher than the r-band CAS parameters. 

To do this we carry out a regression analysis of the 
PCA components for each physical property against the 



CAS analysis of the r-band and it-band light distribution 
of the galaxy. The results are shown in Table [I] for the full 
sample and Table [2] for each of the smaller redshift range 
subsamples. The values in Table [T] and Table [2] show the cor- 
relation coefficients between the distribution of each physical 
property in a galaxy (described by 8 principal components) 
and the three CAS parameters. They also show the correla- 
tion coefficients between the r-band and u-band light distri- 
bution (also described by 8 principal components) and the 
CAS parameters. 

For this analysis, the CAS parameterization of the r- 
band images provides a convenient proxy for a single-filter 
derived morphology. Comparing the QMM results to the 
r-band CAS parameters, as well as to the Pixel-z : CAS out- 
put (physical properties), allows us to establish whether 
the QMM measurements contain more information about 
the distribution of physical properties within a galaxy than 
about the simple distribution of light. 

If the correlation of the PCA components is higher for 
the CAS analysis of the physical properties compared to 
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Figure 4. Principal component analysis for galaxies in the redshift bins. The distributions are very similar and a statistical analysis of 
these similarities is shown in Table 2. 



Table 1. Correlation coefficients between the CAS parameters and the distribution of physical properties in the galaxy (including the 
r-band and u-band light distributions. 



Physical Properties Light Profiles 
(jAo 1 <ir<uiie1 ers 



Age SFR Colour Excess Mctallicity r-band u-band 



Concentration 0.04 0.01 0.02 0.03 0.20 0.04 

Asymmetry 0.18 0.24 0.21 0.20 0.27 0.21 

Clumpiness 0.30 0.46 0.43 0.46 0.28 0.40 



the r-band light profile, QMM would seem to be better at 
representing the spatial distribution of physical properties 
in galaxies than just their simple r-band light distribution. 
The regression analysis was carried out on the entire galaxy 
sample as well as independently for the redshift subsamples, 
in order to identify any systematic effect due to varying red- 
shifts. The tables report the correlation coefficients between 



the PCA output and each of the CAS parameters for each 
property, and for the r-band light distribution. 

There is essentially no correlation between the concen- 
tration of the physical properties and the principal compo- 
nents, particularly when these values are compared with the 
correlation for the r-band light distribution. The asymme- 
try of physical parameters has a slightly higher correlation 
with the output of QMM but even this is weaker than the 
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correlation between QMM and the asymmetry of the r-band 
light distribution. Thus, the QMM morphology cannot de- 
scribe the concentration and asymmetry of physical parame- 
ters better than it can describe the r-band light distribution. 

The correlations between the physical properties in the 
dumpiness parameter and the QMM results are much higher 
than the correlation between QMM and the r-band light dis- 
tribution. This indicates that the results of QMM provide 
a more accurate representation of the underlying physical 
parameter distribution, for the dumpiness parameter, than 
of the r-band light distribution of the galaxy. Thus QMM, 
a morphology indicator, can represent the dumpiness mea- 
surement of physical properties of galaxies better than the 
single-filter photometric light distribution. 

However, there is a possibility that the higher correla- 
tion we find between the dumpiness parameter and QMM 
output could be the result of noise in the tt-band data from 
the SDSS images which was also processed by the QMM 
technique. This noise would make the image appear clumpy 
hence providing a better correlation with the dumpiness pa- 
rameter. On carrying out the correlation analysis with the 
QMM output and the w-band light distribution it is appar- 
ent that dumpiness of the star formation rate, colour excess 
and metallicity still have a higher correlation with the QMM 
output than the u-band output. Thus we can conclude that 
the correlation between the QMM output and the dumpi- 
ness of physical properties is not due to any noise in the 
it-band images and that QMM actually does represent the 
dumpiness of physical properties. 

Following the analysis for the complete sample, the sub- 
samples in redshift bins were analyzed independently, with 
similar results, confirming that there is little or no system- 
atic effect based on redshift (see Table [2b. 

For the separate redshift bins, all show similar patterns, 
and all are in agreement with the results for the complete 
sample. These results confirm that QMM is a better descrip- 
tor of the dumpiness of the physical properties compared to 
the r-band light distribution of the galaxies. For the concen- 
tration and asymmetry parameters QMM cannot describe 
the distribution of the physical properties any better than 
it describes the r-band light distribution of the galaxies. 

An important finding here is that the relative correla- 
tions between physical properties and the r-band light distri- 
butions are consistent across all redshift ranges. The actual 
values between different redshifts vary, but they are still con- 
sistent with the findings from the complete sample. These 
fluctuations can be attributed to the fact that the redshift 
ranges themselves are not even. The two middle bins are 
particularly small compared to the two redshift ranges at 
the edges of the sample. The ranges were chosen such that 
the number of galaxies in each is equal, so as to ensure that 
the principal component analysis is not biased by different 
sample sizes. Furthermore, the fluctuations do not display 
a trend with higher redshifts, indicating that there is no 
real systematic effect due to redshift. Thus we conclude that 
QMM seems to provide a useful approach toward connect- 
ing galaxy morphologies to the morphologies of the physical 
parameters, and is independent of systematic effects due to 
redshift, over the small range considered here. 

We stress that the correlation coefficients should not be 
taken as indicating real correlations. The correlation coeffi- 
cients in all cases are too small to indicate any strong corre- 



lations between the tested parameters. Our goal is to com- 
pare the strength of the coefficients in relation to the r and 
it-bands and identify those CAS parameters for the physical 
properties that have consistently higher correlations, if any, 
than in the r and u bands. 



6 DISCUSSION 

We have two main findings. The first is that the correlation 
between the results of QMM and the dumpiness of physical 
properties of Pixel-z : CAS is significantly higher than the 
correlation with the dumpiness of r-band or u-band light 
distributions. For the concentration and asymmetry param- 
eters, the correlation between QMM and the physical prop- 
erties is much lower than the correlation with the r-band 
light distribution. 

It is expected that the QMM results will be better cor- 
related with the r-band light distribution of the galaxy as 
compared to the derived physical quantities as the QMM re- 
sults are calculated from the light distribution of the galaxy. 
Thus it is indeed significant that the QMM results are bet- 
ter correlated with the dumpiness of the derived physical 
properties than with the dumpiness of the r-band light dis- 
tribution. 

The lack of correlation in the concentration and asym- 
metry parameters is perhaps not completely unexpected. 
Neither of these parameters are particularly effective at de- 
scribing the "patchy" nature of the distribution of physi- 
cal properties in galaxies. The Pixel-z decomposition of the 
spatial distribution of various physical properties in galax- 



ies (as seen in Welikala et al. 2008 20091, is an excellent 



example of this type of distribution. Unlike the light profile 
of galaxies, physical properties are not necessarily centrally 
concentrated. As a consequence, it seems the concentration 
parameter is unlikely to be a useful descriptor for charac- 
terising the distribution of physical properties of galaxies, 
particularly compared to the optical light distribution. 

The asymmetry parameter may also not be effective at 
describing the "patchy" distribution of physical properties, 
and it is further limited given that a center of rotation needs 
to be identified. The brightest pixel was used as the center 
of rotation, which is a reasonable approximation for the op- 
tical light distribution, but the light and physical properties 
do not share well correlated distributions and the high spa- 
tial frequency nature of the physical properties means that 
asymmetry is probably not the best approach for quantifying 
these. Moreover, such "patchy" distributions are inherently 
asymmetric, and many systems would show an asymmetry 
close to the maximum possible asymmetry value (Asymme- 
try^). For these reasons asymmetry, too, seems not to be a 
particularly useful parameter in describing the distribution 
of physical properties in galaxies. 

Finally, there is the dumpiness parameter, which would 
clearly seem to be a more useful statistic in quantifying these 
"patchy" distributions of physical properties of galaxies than 
the previous two, as it was developed to measure this type 
of distribution. The correlation coefficients measured do not 
indicate a real correlation, but we serve to distinguish the 
relative correlation strength of QMM with the pixel-z : CAS 
measurements compared to r and u-bands. The range of 
values of dumpiness is more likely to be encoding much 
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Table 2. CAS parameter values for the physical properties across each redshift bin as well as the r-band light distribution. 



Rcdshift Range CAS Parameters 



Physical Properties 



Age SFR Colour Excess Mctallicity 



Light Profiles 



r-band u-band 



0.00278s; z <0.066 



Concentration 
Asymmetry 
dumpiness 



0.05 
0.20 
0.26 



0.16 
0.27 
0.46 



0.04 
0.23 
0.42 



0.04 
0.23 
0.44 



0.23 0.05 
0.30 0.21 
0.24 0.37 



0.067s; z <0.075 



Concentration 
Asymmetry 
dumpiness 



0.04 
0.20 
0.31 



0.02 
0.27 
0.46 



0.02 
0.21 
0.42 



0.03 
0.20 
0.46 



0.20 
0.28 
0.28 



0.05 
0.21 
0.40 



Concentration 0.22 0.02 0.20 0.03 0.19 0.06 

0.076s: z <0.083 Asymmetry 0.18 0.24 0.20 0.18 0.26 0.21 

dumpiness 0.32 0.47 0.43 0.47 0.36 0.39 



Concentration 0.19 0.11 0.02 0.26 0.19 0.06 

0.084s; z <0.231 Asymmetry 0.17 0.24 0.20 0.19 0.26 0.28 

dumpiness 0.32 0.45 0.42 0.46 0.34 0.49 



more physical information about galaxies, rather than the 
less useful and potentially biased estimates for concentra- 
tion and asymmetry. The QMM morphology should also be 
sensitive to the "patchy" distributions of pixel colours that 
underlies the clumpy nature of the physical properties. This 
could well be the reason for the higher correlation values 
with QMM seen for dumpiness compared to the concentra- 
tion and asymmetry as well as the dumpiness of the r-band 
light distribution. We also confirmed that this dumpiness 
is not due to any noise in the u-band images by showing 
that the correlation between the dumpiness parameters and 
QMM output is higher than that between the it-band and 
the QMM output. This was true for all redshift bins as well. 
Using other approaches towards quantifying clumpy struc- 
ture may be helpful in confirming this. 



7 SUMMARY 

We have compared two separate techniques for analysing 
galaxies, one using a purely physical approach (Pixel- 
z : CAS) and the other using a purely morphological ap- 
proach (QMM). The purely physical approach uses the 
Pixel-z method to infer the distribution of four physical 
properties (age, star formation rate, colour excess and metal- 
licity), within galaxies. These distributions were then quan- 
tified using the CAS method where each physical prop- 
erty was assigned a value for its concentration, asymme- 
try and dumpiness. This is the first application of CAS to 
distributed physical properties within galaxies. The Pixel- 
z : CAS procedure was also applied to the r-band and u-band 
light distribution. We did not discover a trivial relationship 
between the optical light distribution of galaxies and the 
distribution of the physical properties in CAS space. This 
places further significance on our objective of using QMM to 
connect the multiwavelength photometric morphology with 



the underlying physical properties. This is not unexpected 
as CAS is able to be applied only to single-filter imaging, 
rather than a full multiwavelength morphology. 

We also analyzed the morphology of the galaxies using 
QMM where the images of the galaxies in five filters were 
decomposed using shapelets, followed by a Principal Com- 
ponent Analysis. 

To measure the possible correlations between Pixel- 
z : CAS and QMM, we carried out a regression analysis for 
the CAS parameters of each physical property against the 8 
principal components and compared these to a similar anal- 
ysis using the CAS parameters for the r-band and u-band 
light distributions. The regression analysis shows the extent 
to which the spatial distribution of physical parameters and 
QMM are related. 

For the concentration and asymmetry parameters no 
correlation was seen. Thus QMM cannot describe the distri- 
bution of physical properties any better than it can describe 
the r-band light distribution for these two parameters. 

The dumpiness of physical parameters is clearly con- 
nected better to QMM than is the dumpiness of the r-band 
or u-band light distributions. This shows that QMM (which 
was developed to quantify morphology) can describe at least 
some aspects of the distribution of physical properties better 
than the single filter morphology of a galaxy. This demon- 
strates the fact that QMM incorporates physical properties 
as well as the morphology of galaxies through its inclusion 
of colour information. 

Variation in redshift did not significantly affect the fi- 
nal results. The relative correlation between the QMM and 
physical properties, and QMM and the optical light distribu- 
tion is consistent across all the redshift bins. The correlation 
between the QMM and physical properties, and QMM and 
the u-band light distribution was also consistent across all 
the redshift bins. 

This is a significant step in the process of developing a 
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complete galaxy classification scheme. The QMM approach 
is a useful technique that provides a new way of classifying 
galaxies. It maybe more closely related to the underlying 
physical properties of a galaxy than traditional morphology 
measures. 

This is the first time that there has been an attempt to 
combine the spatial distribution of physical properties and 
morphologies of galaxies. The results are promising, with a 
definite connection between the spatial distribution of phys- 
ical properties and the morphology of the galaxies. This sug- 
gests that further investigation is warranted to explore the 
links between morphology and the underlying physical prop- 
erties of galaxies. 
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